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Abstract. We investigated the kinematics of the pulsar wind nebula (PWN) associated with PSR B1951+32 
in the old supernova remnant CTB 80 using the Fabry-Perot interferometer of the 6 m Special Astrophysical 
Observatory telescope. In addition to the previously known expansion of the system of bright filaments with a 
velocity of 100-200 km s _1 , we detected weak high-velocity features in the Ha line at least up to velocities of 
400-450 km s _1 . We analyzed the morphology of the PWN in the Ha, [SII], and [OIII] lines using HST archival 
data and discuss its nature. The shape of the central filamentary shell, which is determined by the emission in 
the [OIII] line and in the radio continuum, is shown to be consistent with the bow-shock model for a significant 
(about 60°) inclination of the pulsar’s velocity vector to the plane of the sky. In this case, the space velocity of 
the pulsar is twice higher than its tangential velocity, i.e., it reaches ~ 500 km s _1 , and PSR B1951+32 is the 
first pulsar whose line-of-sight velocity (of about 400 km/s) has been estimated from the PWN observations. The 
shcll-like Ha-structures outside the bow shock front in the east and the west may be associated with both the 
pulsar’s jets and the pulsar-wind breakthrough due to the layered structure of the extended CTB 80 shell. 
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INTRODUCTION 

CTB 80 is a classical example of a supernova remnant 
(SNR) with a fast-moving pulsar at a late stage of the 
pulsar’s interaction with a very old shell. 

The radio image of CTB 80 is represented by three ex¬ 
tended (about 30') ridges that converge in the region of a 
bright compact core (Velusami and Kundu 1974; Velusami 
etal. 1976; Angerhofer etal. 1981; Strom etal. 1984; 
Mantovani etal. 1985; Strom 1987; Castelleti etal. 2003, 
and references therein). This unusual (for SNRs) morphol¬ 
ogy ceased to appear puzzling in 1988, when, on the one 
hand, the 39.5-ms pulsar PSR B1951+32 was discovered in 
the core (Kulkarni et al. 1988; Fruchter et al. 1988) and, on 
the other hand, Fesen et al. (1988) identified an extended 
infrared shell that extends the radio ridges in the northeast 
and assumed that precisely this infrared shell represents 
a very old SNR. The radio ridges with a steep spectrum 
(a ~ —0.7) correspond to the part of the shell into the 
compressed magnetic field of which the approached pul¬ 
sar injected fresh relativistic particles, thereby reanimat¬ 
ing its synchrotron radio emission. Subsequently, an HI 
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shell that coincides with the infrared shell and expands 
with a velocity of 72 km s _1 , which yields a kinematic 
SNR age of 7.7x 10 4 yr, was also identified (Koo et al. 1990, 
1993). The characteristic age of the pulsar PSR B1951+32 
is t ~ 10 5 yr (Kulkarni etal. 1988; Fruchter etal. 1988). 

The pulsar wind nebula (PWN) is represented by a 
bright compact (45") core with a flat spectrum (a ~ 0.0). 
The core lies at the western boundary of a 10' x 6' plateau 
with a steeper spectrum (a ~ —0.3) elongated in the east- 
west direction. Recently, Migliazzo et al. (2002) detected 
the motion of the pulsar with a velocity of 240 km s~ 4 in 
a direction that confirms the possibility of its birth inside 
the infrared and HI shells. 

The stage of pulsar wind interaction with a very old 
shell observed in CTB 80 seems most complex in the evo¬ 
lutionary chain of PWNs, from young SNRs where the 
pulsar interacts with the SN ejection to pulsars that have 
already escaped from the old SNR. The difficulty lies in the 
fact that the matter density in the cooled old shell behind 
the front of the decelerated shock (and, hence, the density 
of the compressed interstellar magnetic field) is several 
hundred times higher than the ambient density and the 
shell structure is unpredictable in advance. In addition, 
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the rotational energy loss by the pulsar PSR B1951+32 is 
large (E = 3.7 x 10 36 erg s _1 , as estimated by Kulkarni 
etal. 1988), which is comparable to the energy input from 
the young Vela X pulsar. 

Various distance estimates for CTB 80 and the pul¬ 
sar PSR B1951+32 lie within the range 1.5-2.5 kpc (see 
Koo et al. 1993; Strom and Stappers 2000; and references 
therein); we use the universally accepted distance of 2 kpc. 

In this paper, we present our observations of the PWN 
in CTB 80 with the 6 m Special Astrophysical Observatory 
(SAO) telescope and analyze narrow-band observational 
data from the HST archive. 

OBSERVATIONS AND DATA REDUCTION 

Interferometric Observations on the 6 m Telescope 

The core of CTB 80 was observed with the 6 m tele¬ 
scope as part of a program entitled “The Kinematics of 
Matter in Pulsar Wind Nebulae” (the main applicant 
is Yu.A. Shibanov). Our interferometric observations of 
CTB 80 were performed on October 9-10, 2001, at the 
prime focus of the 6 m telescope using the SCORPIO 
focal reducer; the equivalent focal ratio of the system was 
F/ 2.9. A description of SCORPIO is given in the paper 
by Afanasiev and Moiseev (2005) and on the Internet 
(http: / / www.sao.ru/hq/moisav/scorpio/scorpio.html); 
the SCOPRIO capabilities in interferometric observations 
were described by Moiseev (2002). The seeing during the 
observations varied within the range l , .'2-2 , . / 2. We used a 
scanning Fabry-Perot interferometer (FPI) in the 235th 
order at the wavelength of the Ha line; the separation be¬ 
tween the neighboring orders of interference, AA = 28 A, 
corresponded to a region free from order overlapping 
of 1270 km s _1 on the radial velocity scale. The width 
of the instrumental FPI profile was FWHM « 2.5 A 
or ~ 110-120 km s . Premonochromatization was 

performed using an interference filter with a half-width of 
AA = 14 A centered on the Ha line. The detector was a 
TK1024 1024 x 1024-pixel CCD array. The observations 
were carried out with 2 x 2-pixel binning to reduce the 
readout time. In each spectral channel, we obtained 
512 x 512-pixel images; at a 0.55"/pixel scale, the total 
field of view was 417. We obtained a total of 32 interfero- 
grams at various FPI plate spacings, so the width of the 
spectral channel was <5A = 0.87 A or 40 km s _1 near Ha. 
The exposure time was 240 s per channel. 

We reduced the observations using software running 
in the IDL environment (Moiseev 2002). After the pri¬ 
mary reduction (debiasing and flat fielding), the obser¬ 
vational data were represented as 512 x 512 x 32-pixel 
data cubes; here, a 32-cliannel spectrum corresponds to 
each pixel. Major difficulties in the data reduction process 
arose when night-sky emission lines were subtracted from 
interferograms. All bright image features were masked, 
and an azimuthally averaged radial sky emission line pro¬ 
file was constructed from the remaining areas; this profile 
was subtracted from the corresponding frame in the data 


cube. This technique allows the sky line intensity varia¬ 
tions during two-hour observations to be effectively cor¬ 
rected. However, almost the entire field of view near the 
core of CTB 80 proved to be filled with weak emission fea¬ 
tures, which was also confirmed by our deep Ha images of 
the PWN. Therefore, the problem of choosing a sufficient 
number of “clean” sky areas in a region of almost 5' in size 
arose. As a result, the weakest features in the field of view 
could be “oversubtracted”, and the Ha emission line pro¬ 
file is severely distorted at these locations. Our estimates 
indicate that the line distortion due to the background 
subtraction may be disregarded for regions in which the 
intensity of the Ha line exceeds 150 photoelectrons per 
CCD pixel, which corresponds to a surface brightness of 
8.6 x 10~ 17 erg s _1 /D (uncorrected for the interstellar 
reddening). 

The large number of background stars in the field al¬ 
lowed us to measure and correct the atmospheric trans¬ 
parency and seeing variations in each image. The resulting 
seeing was 2'.'2. Subsequently, we reduced the spectra in 
the data cube to the same wavelength scale. The formal 
accuracy of measuring the relative line-of-sight velocities 
was about 2-4 km s _1 . However, since the emission lines 
in the object often have asymmetric profiles, the actual 
measurement accuracy depends on the chosen line profile 
fitting method. 

The spectra of the object were smoothed with a 
FWHM = 1.7 A (two channels) Gaussian for optimal 
filtering of the data cube. To reliably identify weak emis¬ 
sion features, we also smoothed the images in the cube 
with a bivariate Gaussian. The resultant angular resolu¬ 
tion was about 2 / /6. The smoothing was performed using 
the ADHOC software package 1 * * . 

The continuum level in each spectrum was determined 
as a median mean of the eight weakest levels. When con¬ 
structing the velocity field and the monochromatic im¬ 
age in the Ha line, we fitted the spectral line profile by a 
Gaussian using only the points offset by ±3 channels from 
the predetermined line peak. 

We calibrated the monochromatic image in energy 
units (erg s^ 1 cm -2 ) by comparison with our calibrated 
narrow-band image (see below). 

All of the radial velocities in this paper are heliocentric; 
the passage to the Local Standard of Rest corresponds to 

y Hel = V LSR + 17 - 6 km s_1 - 

Medium- and Narrow-Band Images 

Observations with the 6 m telescope. Our mediurn- 
and narrow-band observations of the CTB 80 core were 
performed on October 11 and 12, 2001, with the 6 m tele¬ 
scope using the SCORPIO focal reducer (see above) at 
seeing l" 4 1"6 at zenith distance 2 = 12°-21°. Filter pa¬ 
rameters and exposure times are given in the Table. 

1 The ADHOC software package was developed by 

J. Boulestex (Marseilles Observatory) and is freely accessible 

on the Internet. 
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Details of the narrow- and medium-band observations of the CTB 80 core with the 6 m telescope (BTA) and HST 


Range, telescope 

Filter 

Acem A 

FWHM, A 

Exposure time, s 

Ha, BTA 

FN657 

6578 

75 

4 x 300 

Continuum, BTA 

SED707 

7036 

207 

4 x 120 

Ha, HST 

F656N 

6563 

20 

2600 + 2700 

Continuum, HST 

F547M 

5479 

200 

1300 + 1300 

[OIII], HST 

F502N 

5013 

20 

2700 + 2700 

[SII], HST 

F673N 

6732 

20 

2700 + 2700 


To calibrate the images, we observed the spectropho- 
tometric standards G138-31 and Feige 110. 

The HST archive. The PWN in CTB 80 was ob¬ 
served with HST (the main applicant is J. Trauger) in 
October 1997 using the WFPC2 instrument. Parameters 
of the filters used and the total exposure time are given 
in the Table. The primary data reduction is automat¬ 
ically performed when the data are queried from the 
HST archive. A filtering code proposed by N.A. Tikhonov 
(private communication) was used to remove numerous 
cosmic-ray particle hits. The formal accuracy of the as¬ 
trometric referencing using WCSTools is an order of mag¬ 
nitude lower than the internal accuracy of the USNO-A2.0 
star catalog used as a reference one. The resultant value 
was taken to be 073. 

RESULTS OF OBSERVATIONS 

The PWN Morphology 

Previous studies of the CTB 80 core showed that its op¬ 
tical emission is typical of SNRs: an intense (relative to 
Ha) [Nil], [SII], [01], [OIII] line emission and a filamen¬ 
tary structure (see Angerhofer et al. 1980; Blair et al. 1984; 
Whitehead et al. 1989; Hester and Kulkarni 1988, 1989; 
and references therein). Significant differences in the PWN 
morphology in different lines were reported: only a sym¬ 
metric central shell is observed in the [OIII] line, a shell 
to the east of it appears in the [SII] line, and two bright 
shells to the east and the west of the central structure that 
form an elongated core structure are observed in the Ha 
line. 

The deep HST images of the PWN demonstrate a stag¬ 
gering filamentary, irregularly shaped multishell structure 
and confirm the results of ground-based observations (see 
Fig. 1, which shows the [OIII], [SII], and Ha images, from 
top to bottom panels, respectively, with superimposed 
VLA 1.5-GHz radio isophotes). 

Only the central horseshoe-shaped part of the core 
closest to the pulsar is seen in the [OIII] line. The PWN 
morphology in the [OIII] line clearly shows a structure 
expected for the bow shock produced by the pulsar’s mo¬ 
tion (see the discussion). The Ha emission is also ob¬ 
served in this central filamentary horseshoe-shaped shell; 
the Ha filaments are adjacent to the [OIII] filaments from 


the outside. We clearly see the two bright (in Ha) fila¬ 
mentary shells in the east and the west, which form a 
spindle-shaped core structure in the east-west direction. 
The PWN sizes in the Ha line are about 75" x 38" or 
0.73 x 0.37 pc. 

The PWN Kinematics 

The velocities of PWN filaments. The comparisons 
of relative line intensities in the core spectrum with di¬ 
agnostic models made by several authors are suggestive 
of collisional excitation of the gas behind the front of a 
shock propagating at a velocity of about 120-140 km s -1 
in a medium with a density of 25-100 cm -3 (Hester and 
Kulkarni 1989, and references therein). Gas motions in 
the PWN with such velocities have already been detected. 
The main method used was long-slit spectroscopy, which 
allowed the velocities of the individual, generally bright¬ 
est (and, hence, possibly slowest) knots and filaments 
to be estimated (Angerhofer etal. 1981; Blair etal. 1984). 
Spectroscopic methods did not reveal radial velocities 
higher than 200 km s _1 with a surface brightness at 
Ha > 0.5 10 6 erg cm 2 s 1 sr 1 anywhere in the PWN 
(Blair etal. 1988). Using the echelle spectrograph of the 
4-m WHT telescope, Whitehead etal. (1989) found an ex¬ 
pansion of the two central shells with velocities of 200 and 
100 km s _1 . 

Using FPI, we have studied the PWN kinematics for 
the first time. The advantage of FPI observations is that 
they give the line profile everywhere in the PWN, and not 
only in the region cut out by the spectrograph slit. 

Our observations indicate that the line profile, which 
is single in bright peripheral filaments, is characterized 
by a multicomponent structure or line asymmetry, bright 
red wings in the central PWN regions. As an illustration, 
several profiles and their decomposition into individual 
Gaussians are given in Fig. 2. Having analyzed all of the 
observed Ha profiles, we were able to localize the high- 
velocity gas in the PWN image that emits in the range 
from 200 to 400 km s _1 (shown at the center in Fig. 2). 

Figure 3 shows several so-called position-velocity di¬ 
agrams (the change in gas velocity along the chosen di¬ 
rection) constructed from our FPI observations and the 
localization of the corresponding directions in the PWN 
image. 
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Fig. 1. HST images of the PWN with superimposed VLA 
1.5-GHz radio isophotes: in the [OIII] line (top), in the 
[SII] line (middle), and in the Ha line (bottom panel). The 
arrow indicates the motion of the pulsar over a period of 
1000 yr, as measured by Migliazzo etal. (2002). 


As we see from this figure, in the bright core structures, 
the position-velocity diagrams show velocities that agree 
with the previously measured expansion velocity of 100- 
200 km s -1 ; i.e., they completely confirm the results of 
spectroscopic observations. 

In addition to the gas velocities in bright filaments 
within the range from —200 to +200 km s^ 1 , our obser¬ 
vations clearly reveal weaker high-velocity Ha emission 
features in the PWN. Weak emission at high velocities is 
clearly seen in the position-velocity diagrams up to 400- 
450 km s -1 . This is a lower limit, since the velocity range is 
limited by the FPI velocity range free from order overlap¬ 
ping and the impossibility to properly take into account 
the contribution from the [Nil] 6548 A line emission. 

Figure 3 (scan 3) clearly shows the typical (for an ex¬ 
panding shell) structure of the so-called velocity ellipsoid 
(its half corresponding to positive velocities) that is deter¬ 
mined by high-velocity motions. The corresponding pos¬ 
sible expansion velocity reaches the lower limit mentioned 
above. 

The surface brightness of the detected high-velocity 
component of the Ha line in the central part of the core lies 
within the range (1.4-20) x 10 -16 erg s -1 cm -2 arcsec~ 2 . 
In our estimation, we took the color excess E(B—V) = 0.8 
from Blair etal. (1984). 

Figure 4 shows the line-of-sight velocity held corre¬ 
sponding to the peak of the main line component con¬ 
structed from our observations and superimposed on the 
Ha image of the PWN. We emphasize that the line-of- 
sight velocity held refers only to the core of the line, which 
often has a complex multicomponent prohle. The veloci¬ 
ties determined from the line peak lie within the range 
— 100 to +50 km s -1 , in agreement with the spectroscopic 
observations. 

In Fig. 4, we clearly see the symmetry axis in the veloc¬ 
ity distribution of the peak of the main line component in a 
direction P ss 230°, which is close to but does not coincide 
with the direction of the pulsar’s motion (P = 252° ± 7°; 
Migliazzo etal. 2002). 

We estimated the total hux from the PWN in the main 
line component to be 8.4 x 10 -12 erg s _1 cm -2 , in close 
agreement with the estimate of Whitehead etal. (1989). 
The total flux in the high-velocity line component is 6.4 x 
10 -13 erg s _1 cm -2 ; the luminosity of the high-velocity 
gas accounts for about 7% of the total luminosity of the 
core in the Ha line. 

The high-velocity emission detected in the core for the 
first time has confirmed the changes in the PWN fine 
structure pointed out by Strom and Blair (1985) — pos¬ 
sibly the proper motions of the filaments of the central 
shell corresponding to a velocity of 250d km s _1 (up to 
400d km s _1 ), where d is the distance to CTB80 in kpc, 
i.e., about 500 km s -1 (possibly up to 800 km s _1 ). 

The velocities of filaments in the CTB 80 shell 
outside the core. The optical radiation from the ex¬ 
tended SNR CTB 80 outside the PWN is represented by 
the system of faint filaments clearly seen in the [SII] and 
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Fig. 2. Ha profiles at several points of the PWN and their Gaussian decomposition. The features at a velocity of 
+605 km s _1 are attributable to the [Nil] 6548 A line emission that falls in the wings of the interference filter and 
that is separated from the Ha emission by about 0.5 order of interference. The localization of the high-velocity gas 
emitting in the range from 200 to 400 km s -1 is shown in the central Ha image of the PWN (isophotes). 


Ha lines over the entire 16' x 16' field whose image is given 
in the paper by Hester and Kulkarni (1989). 

The filamentary morphology and the intense [SII] line 
emission provide strong evidence for the gas emission be¬ 
hind the shock front. It is interesting to measure the ve¬ 
locities of these filaments in the extended remnant out¬ 
side the core, since, in general, the thin filaments can be 
localized on the front or rear side of the old shell; i.e., 
these may be projected rather than be physically associ¬ 
ated with the PWN. Besides, the pulsar’s velocity with 
respect to the ambient gas is important in estimating the 
expected shape of the bow shock produced by the wind 


from a fast-moving pulsar. Figure 4 (bottom panel) shows 
only the brightest outer filaments in a 4' x 4' field for 
which the background subtraction effect is insignificant 
(see above). The Ha brightness of these filaments exceeds 
8.6 x 1CU 1 ' erg s^/D . As follows from Fig. 4 (bottom), 
the velocities of the bright filaments in the outer shell of 
CTB 80 do not differ significantly from the velocities of 
the peripheral filaments of the PWN undistorted by the 
expansion of the latter. 

Note, in particular, the filaments that are immediately 
adjacent to the PWN in the north and the east and that 
form a kind of a wake at the boundary of the shell pro- 
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Fig. 3. Position-velocity diagrams constructed from our FPI observations and the localization of the corresponding 
scans in the PWN image. The emission features at a velocity of about 600 km s -1 are attributable to the [Nil] 6548 A 
line. 


duced by the wind from a moving pulsar. (These filaments 
are most distinct in the [SII] lines in Fig. 2b from Hester 
and Kulkarni 1989). As our measurements indicate, the ve¬ 
locities of these filaments do not differ significantly from 
those of the bright PWN filaments either. 


DISCUSSION 

The difficulty of explaining the nature of the PWN in 
CTB 80 lies in the fact the bow-shock structure (in the 
[OIII] line, the radio continuum, and the soft X-ray band) 
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Fig. 4. CTB80 core. Top panel: the velocity field of the 
Ha peak superimposed on the PWN image. The arrow 
indicates the direction of the pulsar’s motion, as measured 
by Migliazzo etal. (2002); bottom panel: the same for the 
entire FOY of about 5' x 5' in size in the image obtained 
with the 6 m telescope (only the brightest outer filaments 
are shown). 


and the elongated filamentary shell-like structure in the 
Ha line, which is definitely also of an shock origin that, in 
the opinion of Hester (2003), is associated with the action 
of the pulsar’s jets, are simultaneously clearly seen around 
the pulsar. We do not see any manifestations of the inter¬ 


action between two processes: the ejections of relativistic 
pulsar plasma in two jets and the motion of the pulsar. 
The elongated PWN structure is not distorted by the pul¬ 
sar’s motion in any way; the bow shock, in turn, is not 
distorted by the action of the two plasma jets in any way 
either. The rough estimates given below indicate that in¬ 
dividually both structures can be explained adequately. 


Estimates for the Shock Waves Produced by the 
Pulsar’s Jets 


If the pulsar’s power L is constant over time t and is re¬ 
leased in a solid angle H, blowing a bubble of radius R = 
R(t) in a homogeneous interstellar medium with a den¬ 
sity po, then the mass M ~ flp 0 R 3 /3 is swept up from 
the bubble. Assuming that all of this mass is gathered in 
a thin layer at radius R and remains in the same solid 
angle H (in fact, it will partially flow around the jet) and 
equating the pulsar-wind momentum flux L/c to the rate 
of increase in the momentum of the wind-driven shell, 


L 

c 


j t {Mv) 


where v = dR/dt , we obtain an acceptable size of the 
bubble swept up by the wind for a pulsar at rest under 
reasonable assumptions about the ambient density and the 
wind asphericity. Assuming that R oc t a , we have 

l n d {r 4 \ 

7 _ 3 ' p ° a Jt \T) ’ 


hence, a = 1/2 at L = const, i.e., 


R = 



1/4 

f l/2 


3 pc x ( £L. 


J 36 


47T 
Hno ) 


1/4 


,1/2 


(1) 


where £ = L/Lq is the fraction of the power Lq went into 
the solid angle Q and spent on blowing the bubble, no is 
the particle number density (for hydrogen composition), 
and tg is the time in Myr. 

The power of the flux in relativistic particles (including 
all types of photons) from the pulsar PSR B1951+32 is 
Lq = -E/ot = 3.7 x 10 36 erg s _1 (Kulkarni etal. 1988). 
At no = 1 cm -3 , L 36 = 3.7, and the pulsar’s total age 
t 6 = 0.1, we obtain a bubble size of about 1 pc at £ = 1 
even for isotropic radiation (H = 47t), in agreement with 
the observed PWN size in the Ha line. In fact, only a 
small fraction £ of the pulsar’s power goes into blowing 
the bubble, but H is also much less than 47r (by several 
hundred times for a linear jet opening angle of several 
degrees). 

We emphasize that this is the minimum upper limit: 
here, the shock waves are assumed to be radiative. The 
bubble size will increase if the hot-gas pressure inside the 
bubble is taken into account (Castelleti etal. 2003). Two 
cases are possible: the radius for radiative shock waves is 
larger than our estimate, but smaller than that for adia¬ 
batic shock waves (McKee and Ostriker 1977; Blinnikov 
etal. 1982). 
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If we took into account the hot-gas pressure inside the 
bubble and for adiabatic shock waves, we would obtain the 
maximum estimate at 0 = 47 t (Avedisova 1971; Weaver 
etal. 1977; Ostriker and McKee 1988): 

R = 28 pcxn 0 - 1/5 L^ /5 . (2) 


To estimate the shock shape, we can use the solutions by 
Lipunov and Prokhorov (1984) and Wilkin (1996). The 
latter author obtained an analytical solution for a normal 
(nonrelativistic) stellar wind with radiative shock waves: 

i?(0) = c?oCSc[3(l — ^cotfi 1 )] 1 ^ 2 , (3) 


Hence, R could be even an order of magnitude larger than 
that from (1). 

Let us make several remarks suggesting that these 
rough estimates are uncertain. On the one hand, the 
B aimer line emission in PWNs is evidence of nonradia- 
tive shock waves in a partially neutral medium (Chevalier 
and Raymond 1980); it is produced by electron impact 
excitation or ion charge exchange. In this case, neutral 
hydrogen can penetrate into the hot gas and impart odd 
shapes to the shock waves (Bucciantini and Bandiera 2001; 
Bucciantini 2002; D’Amico etal. 2003) 

On the other hand, formula (2) was derived for a hot 
wind with a normal (Pascal) pressure. Since the flow of 
magnetized particles or photons in a pulsar wind is di¬ 
rected only along the radius, a simple formula of form (1) 
the derivation of which assumes no pressure isotropy holds 
good, while formula (2) overestimates the result. Many 
problems associated with the composition of this wind 
have not yet been solved (the so-called sigma paradox 
— the transition from the fraction of the electromagnetic 
pressure expressed in terms of the Poynting vector to the 
kinetic pressure of the particle flux cannot be reliably cal¬ 
culated (see the review by D’Amico etal. 2003)). 

A time shorter than the pulsar’s age should be taken 
as te , since it has flown into the dense layers of the shell of 
the old SNR relatively recently. Castelleti et al. (2003) took 
t = 18 200yr as the PWN age at an ambient density of 
no = 0.5 cm -3 . According to Mavromatakis etal. (2001), 
the relative line intensities in the spectrum of the filaments 
of the extended CTB 80 shell are typical of the gas radia¬ 
tion behind the front of a shock propagating at a velocity 
of 85-120 km s -1 in a medium with an initial density of 
about 2-5 cm -3 . This density is equal to the mean density 
in the HI shell estimated by Koo etal. (1990). 

However, the dependence of R on all parameters, in¬ 
cluding the density and the age, is weak; i.e., clearly, a 
small fraction t of the total power that goes into blowing 
the bubble will suffice. (We also see the pulsar’s emission; 
i.e., clearly, the fraction of the captured power, £, is ap¬ 
preciably smaller than unity.) 

Thus, the sizes of the bubbles swept up by the pulsar- 
wind jets can easily be made close to the observed values 
even if the density is much higher and the fraction t of the 
pulsar’s power that goes into the jet is small. 

An Estimate for the Location of the Bow Shock 

The pulsar’s motion in the plane of the sky with a velocity 
of Vo = 240 km s -1 (Migliazzo etal. 2002) for an isotropic 
wind must give rise to a bow shock, which is observed 
in several pulsars (see the review by D’Amico etal. 2003). 


where do is the distance from the pulsar to the head point 
at contact discontinuity. This solution describes well the 
classical numerical results by Baranov et al. (1971) and van 
der Swaluw etal. (2003), see Fig.8 in the latter. An ap¬ 
proximate solution for adiabatic shock waves was given 
by Chen etal.{ 1996). The criticism of the applicability 
of such solutions to known PWNs (see Bucciantini and 
Bandiera 2001) is based on the fact that these do not in¬ 
clude the effects of charge exchange and neutral hydrogen 
penetration behind the shock mentioned above. However, 
the PWN in CTB 80 is unique in that the bow shock is 
observed not only in Balmer lines, but also, most clearly, 
in [OIII] lines and in the X-ray and radio bands, i.e., in a 
hot, ionized plasma. Therefore, simple analytical solutions 
may well be valid in this case. 

We obtain the following estimate from the balance be¬ 
tween the momentum fluxes L/c and 47 tg?q/V{;q: 


h L \ 1/2 

d-KCpoVl ) 


( 4 ) 


Hence, we find for vq = 240 km s -1 , p 0 = 2 x 10 -24 g cm -3 
(no = 1 cm -3 ), and £\ = 1 that do ~ 0.044 pc. This value 
is comparable to the observed distance from the pulsar 
to the [OIII] filaments that determine the position of the 
head point at the shock front, d () | )S ~ 6" « 0.057 pc. 

The theoretical value of do will be close to d Q ^ if we 
take into account the fact that no = 170 km s -1 should be 
substituted for vo = 240 km s -1 , since the pulsar moves 
in the matter of a remnant expanding with a velocity of 
72 km s -1 . In this case, do « 0.043 pc. 

In Fig. 5 (top panel), the solution by Wilkin (1996) is 
superimposed on the PWN image in the [OIII] line and 
in the radio continuum. The distance from the head point 
at the front to the pulsar is assumed to be do « 0.043 pc, 
and the axis of the surface with a position angle of 235° 
does not coincide with the pulsar’s velocity direction in 
the plane of the sky. An azimuthal asymmetry in the pro¬ 
jection of the theoretical surface relative to the observed 
shape of the central shell remains noticeable at a posi¬ 
tion angle of P — 252° according to the measurements by 
Migliazzo etal. (2002). 

The fact that the shape of the bow shock observed 
in the [OIII] line is appreciably “broader” than the gas- 
dynamical solutions by Wilkin (1996) may suggest that 
the pulsar moves at a significant angle to the plane of the 
sky. Allowing for the possible inclination of the pulsar’s 
velocity vector to the plane of the sky yields better agree¬ 
ment of the theory with the observed morphology of the 
central shell. By varying the pulsar’s velocity vector, we 
found the best agreement with the observations for a bow 
shock whose axis is inclined at an angle of 60° to the plane 
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Fig. 5. Top panel: PWN in the [OIII] line and in the radio 
continuum (isophotes): the bow shock corresponding to 
solution (3) for the pulsar’s motion in the plane of the sky 
with a position angle of 235° (dashed line) and the bow 
shock inclined at an angle of 60° to the plane of the sky 
(9 isolines). Bottom panel: The same bow shock inclined 
at an angle of 60° chosen by the best agreement with the 
[OIII] line image superimposed onto the Ha and X-ray 
image (isophotes) of the PWN. The arrow indicates the 
pulsar’s motion over a period of 1000 yr, as measured by 
Migliazzo et al. (2002). 


of the sky. Figure 5 (top panel) shows this solution for the 
bow shock chosen by the best agreement with the shape 
of the central shell observed in the [OIII] line and in the 
radio continuum. As we see, the shape of the projection of 
the theoretical surface defined by relation (3) can be com¬ 
pletely reconciled with the observed PWN morphology. In 
this case, the position angle of the pulsar’s velocity vector 
in the plane of the sky is P = 235°; i.e., it is close to the 
direction of the symmetry axis in the velocity field of the 


main Ha- line component found above. P = 235° matches 
the value obtained by Strom (1987). 

We emphasize that, if the pulsar’s velocity vector is 
inclined at an angle of 60°, its space velocity is twice as 
high as its velocity measured in the plane of the sky; i.e., it 
reaches 500 km s _1 (this value agrees with the mean value 
in the velocity distribution of pulsars; see Arzoumanian 
etal. 2002.) Since the distance do at a higher velocity is a 
factor of \fl smaller, for close agreement with the obser¬ 
vations, we must take a slightly lower gas density in front 
of the bow shock. 

Our kinematic studies are also consistent with the sug¬ 
gested model. Since the thickness of the postshock emit¬ 
ting gas, which is determined by the apparent thickness of 
the bright peripheral PWN filaments, is an order of mag¬ 
nitude smaller than the distance the pulsar traverses in 
lOOOyr, (see Fig. 1), we conclude that the radiative cool¬ 
ing time of the postshock gas is short, about 100 yr. In 
this time, an element of gas is not only compressed, but 
also acquires a velocity that corresponds to the pulsar’s 
motion. Without detailed calculations, it is hard to tell 
at which velocity the gas emission is at a maximum; one 
may only expect the observed velocity to be a significant 
fraction of the pulsar’s velocity. 

Gas motions with such velocities in the PWN have 
been detected in our work for the first time (see above). 
These velocities refer to the weak emission features in 
the Ha- line; they are observed both in the central region 
around the pulsar, and near the bright filaments at the 
PWN boundary (see Fig. 2). The line-of-sight velocities of 
the bright filaments are within the ±200 km s -1 range. 
This is an ordinary (for SNRs) situation related to the 
fact that the bright filaments usually represent the front 
surfaces seen edge-on. 

Unfortunately, since the velocity measurement range 
in our FPI observations is limited, at present, we cannot 
unambiguously determine whether the pulsar is moving 
toward or away from us. The existence of high positive 
velocities argues for the motion away from us, but as yet 
we have no information about the possible negative veloc¬ 
ities without further FPI observations. The derived veloc¬ 
ity distribution of the line peak could clarify the picture, 
but here we are restricted by the fact that the density 
distribution in the closest neighborhood of the shell cor¬ 
responding to the bow shock is unknown. In particular, 
the asymmetric expansion of the central shell (according 
to Whitehead etal. (1989), the motion of the approaching 
side at a velocity of about —200 km s -1 is most clearly ob¬ 
served in it) could be attributable to the higher brightness 
of this side due to a nonuniform ambient gas density. 

In Fig. 5 (bottom panel), the same theoretical surface 
of the bow shock that is inclined at an angle of 60° to the 
plane of the sky and that agrees best with the PWN emis¬ 
sion in the [OIII] line is superimposed on the Ha image 
with X-ray isophotes. The two shell-like Ha structures in 
the west and the east are far outside the bow shock. These 
structures, which form the elongated PWN shape, were 
explained by the action of the pulsar’s jets (Hester 2000). 



10 


LOZINSKAYA et al.: NEW STUDIES OF THE PULSAR WIND NEBULA 


It should be noted that, if the pulsar’s jets are directed 
at an angle to its space velocity (e.g., lie in the plane of 
the sky), then the contradiction with the absence of an 
apparent influence of the jets on the shock front in its 
brightest parts mentioned at the beginning of the section 
is removed. 

Of course, the representation of the shock as an in¬ 
finitely thin layer with shape (3) is oversimplified. More 
detailed solutions should be used, and the structure of two 
shocks with a contact discontinuity between them should 
be taken into account (Baranov etal. 1976), since there is 
a hint at such a structure in the observations. In particu¬ 
lar, Moon et al. (2004) associate the structure observed in 
X rays with a backward shock wave in the pulsar wind and 
in the Ha line with a bow shock in the ambient medium. 

Note also that, in fact, i\ in formula (4) may be the 
fraction of the pulsar’s emission absorbed by an ionized 
plasma, while l in formula (1) for the jet may be the 
fraction of the emission absorbed by neutral hydrogen. 
For example, hard ultraviolet emission must be absorbed 
completely by neutral hydrogen and only partially by ion¬ 
ized plasma. Another possible process, the proton-beam 
charge exchange in a neutral gas, cannot take place in an 
ionized medium. Thus, different components of the pul¬ 
sar’s flux can manifest themselves differently in different 
components of the ambient medium. 

A different interpretation of the elongated PWN shape 
in the Ha line is also possible. The system of thin filaments 
mentioned in the Section “Results of Observations”, which 
characterizes the optical emission from the extended rem¬ 
nant outside the PWN, is most likely the layered structure 
of the old CTB 80 shell (Hester and Kulkarni 1989). Our 
velocity measurements for the brightest filaments suggest 
that these are not projected, but are physically associated 
with the PWN. 

The thickness of the layers (typically, ~ 10-20" or 0.1 
0.2 pc) and the separation between them (about 2'-2.5' 
or 1-1.5 pc) at the pulsar’s velocity of 240 km s _1 in the 
plane of the sky yield a time of about 5000 yr between 
the passages of neighboring layers and a passage time of 
a dense layer of 400-800 yr. 

The pulsar’s motion through such a layered medium 
can produce the observed elongated multishell structure 
of the core as a result of the pulsar-wind breakthrough 
from a dense layer into a tenuous medium between the 
layers. The breakthrough of the nonrelativistic wind of 
Wolf- Rayet stars from a dense cloud into a tenuous inter¬ 
cloud gas leads precisely to this effect (see, e.g., Dopita 
and Lozinskaya 1990). 

Localization of the Pulsar in the Extended Remnant 
CTB 80 

If the space velocity of the pulsar actually reaches 500 
km S', then, in general, it traverses a distance of 50 pc 
in a time of 10 5 yr and could go beyond the symmetric 
infrared and HI shells about 40 pc in size, which, accord¬ 


ing to the previous interpretation, determined the total 
volume of the old remnant CTB 80. However, the deep 
images taken in the Ha+[NII], [SII], [Oil], [OIII] lines 
by Mavromatakis et al. (2001) revealed large-scale filamen¬ 
tary and diffuse structures in the nearby 2° x 2° region. 
These optical filaments in the north and the northwest 
closely correlate with the radio images of the remnant 
(Castellcti etal. 2003) and with the boundary of the in¬ 
frared shell, but go far beyond the infrared and HI shells 
in the south and the east. (Here, we disregard the extended 
filaments in the east associated not with CTB 80, but with 
the remnant of a different supernova (Mavromatakis and 
Strom 2002).) The relative line intensities in the spectrum 
of these large-scale filaments are typical of the radiative 
cooling of the gas behind the front of a shock propagating 
at a velocity of 85-120 km s -1 in a medium with an ini¬ 
tial density of about 2-5 cm -3 (Mavromatakis et al. 2001). 
Such a density agrees with the mean density in the HI 
shell estimated by Koo etal. (1990) from 21-crn line obser¬ 
vations. This argues for the association of the large-scale 
filaments with the SNR CTB 80. 

Therefore, a more complex spatial structure composed 
of two hemispheres with different sizes should probably be 
considered as the SNR CTB 80. One hemisphere, which is 
determined by the northeastern radio ridge and the in¬ 
frared and HI shells, is the result of the interaction of 
the shock triggered by a supernova explosion with a dense 
medium. The second part of the shell has an aspherical 
shape and is determined by the southwestern ridge and 
the [OIII] filaments denoted by III and IV as well as by 
the system of bright filaments to the east of IV (see Fig. 2 
from Mavromatakis etal. 2001). This part of the shell was 
most likely formed by the shock in a medium with a much 
lower density. 

Figure 6 shows approximate boundaries of the two 
parts of the shell in projection onto the plane of the sky 
in this model. 

Note that ROSAT observations revealed a conical ~ 1° 
region of thermal X-ray emission southeast of the pul¬ 
sar far beyond the infrared shell (Safi-Harb etal. 1995). 
The central region of the second part of the shell men¬ 
tioned above could, in principle, be responsible for this 
emission. However, the SNR CTB 80 is observed along 
the Cygnus spiral arm and is immediately adjacent to 
the giant Superbubble produced by intense stellar wind 
from the Cyg OB2 cluster (Lozinskaya et al. 2002, and ref¬ 
erences therein). Therefore, in general, the thermal X-ray 
emission in the extended conical region could be the back¬ 
ground emission, i.e., it could belong not to CTB 80, but 
to the Superbubble. 

In the proposed scheme, the morphology of the SNR 
in the plane of the sky suggests that the major axis of 
this structure composed of two hemispheres with different 
sizes is oriented at a large angle to the plane of the sky. 
Therefore, at a possible space velocity of 500 km s _1 , the 
pulsar has not yet gone outside the shell of the extended 
remnant CTB 80. 
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Fig. 6. General scheme of CTB80 from Mavromatakis etal. (2001): the [SII] image of the region, radio ridges (solid 
line), and the infrared and HI shells (the dashed and dash-dotted lines, respectively). The large and small circles 
indicate the suggested model for the shell of the old SNR composed of two hemispheres in projection onto the plane 
of the sky. 


CONCLUSIONS 

Our kinematic study of the pulsar wind nebula in the old 
supernova remnant CTB 80 using the FPI of the 6 m SAO 
telescope has revealed weak high-velocity Ha features in 
the PWN at least up to a velocity of 400-450 km s -1 . We 
confirmed the previously measured expansion of the sys¬ 
tem of bright filaments with a velocity of 100-200 km s -1 . 
We analyzed the PWN morphology in the Ha, [SII], and 
[OIII] lines using the HST archival data. The shape of the 
bow shock, which is determined by the central horseshoe¬ 
shaped shell bright in the [OIII] line and in the radio con¬ 
tinuum, was shown to be in best agreement with the the¬ 
ory for a significant (about 60°) inclination of the pulsar’s 
velocity vector to the plane of the sky. The space velocity 
of the pulsar is twice as high as its tangential velocity mea¬ 
sured by Migliazzo etal. (2002); i.e., it reaches 500 km s . 
This pattern of motion is also confirmed by the high radial 
velocities of the gas in the PWN that we found here. Thus, 


PSR B1951+32 is the first pulsar whose possible light-of- 
sight velocity (of about 400 km s _ l) has been estimated 
from the PWN morphology and kinematics. 

The filamentary shell-like structures observed in the 
Ha line in the east and the west outside the bow shock 
can be explained not only by the action of the pulsar’s 
jets, but also by the pulsar-wind breakthrough into an 
inhomogeneous ambient medium. 

We considered the general scheme of CTB 80 that in¬ 
cludes the most recent optical and radio observational 
data in which the pulsar’s high space velocity is consis¬ 
tent with its location in the dense shell of the old SNR. 

Of course, the proposed scheme of the PWN in CTB 80 
must be confirmed additionally. The existence of high ve¬ 
locities in the PWN, the pattern of elongated Ha struc¬ 
tures (possibly under the influence of the pulsar’s jets or 
when the wind breaks through from a thin dense layer), 
and other questions require further observations and de- 
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tailed hydrodynamic simulations of this interesting object. 
An analysis of the velocity field for the extended filaments 
in the entire region can give a kinematic confirmation for 
the proposed scheme of the old SNR CTB 80. 
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